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ABSTRACT 

We report a discovery of a proto-cluster in vigorous assembly and hosting strong star forming 
activities, associated to a radio galaxy USS 1558-003 at z=2.53, as traced by a wide-field narrow-band 
Ha imaging with MOIRCS on Subaru Telescope. We find 68 Ha emitters with dust-uncorrected SFRs 
down to 8.6 M Q yr _1 . Their spatial distribution indicates that there are three prominent clumps of Ha 
emitters, one surrounding the radio galaxy and another located at ^1.5 Mpc away to the south-west, 
and the other located in between the two. These contiguous three systems are very likely to merge 
together in the near future and may grow to a single more massive cluster at later times. Whilst 
most Ha emitters reside in the "blue cloud" on the color-magnitude diagram, some emitters have 
very red colors with J — K s > 1.38(AB). Interestingly, such red Ha emitters are located towards the 
faint end of the red sequence, and they tend to be located in the high density clumps. We do not see 
any statistically significant difference in the distributions of individual star formation rates or stellar 
masses of the Ha emitters between the dense clumps and the other regions, suggesting that this is one 
of the notable sites where the progenitors of massive galaxies in the present-day clusters were in their 
vigorous formation phase. Finally, we find that Ha emission of the radio galaxy is fairly extended 
spatially over ~ 4.5". However it is not as widespread as its Lya halo, meaning that the Lya emission 
is indeed severely extended by resonant scattering. 

Subject headings: galaxies: clusters: general - galaxies: clusters: individual: USS 1558-003 - galaxies: 
evolution. 



1. INTRODUCTION 



It is well-known that galaxy properties are strongly de- 
pendent on the envi r onment where galaxies reside (e.g.. 
Dressier et ail 119971; iTanaka et ail 120051: iCucciati et ail 



20061 ). Galaxy clusters are one of the most biased en- 
vironments in the Universe, and local ones are dom- 
inated by galaxies with red colors and elliptical mor- 
phologies which have already quenched their star forming 
activities and are passively evolving. In contrast, star 
forming galaxies with a blue color and spiral/irregular 
morphologies are preferentially found in the field envi- 
ronment. This trend is known as a star formation - 
density relation where star formation activity of galax- 
ies decreases gradually in higher density r egions (e.g., 
iKauffmann et all l2004t iCooper et all [2008) . It is cer- 
tain that some processes specific to high density envi- 
ronment play an important role in formation and evo- 
lution of elliptical galaxies found in the local Universe. 
The plausible environmental effects include galaxy- 
intracluster mediu m (ICM) interaction s such as ram- 
pressure stripping (|Gunn fc Gottl 11972ft. gal axy-cluster 
gravitational interactions (jBvrd &; Valtonenlll990T ). and 
galax y-galaxy interactio ns such as mergers and harass- 
ment (M oore et all 119961 ). However, the relative impor- 
tance of the processes still remains unclear. 
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The fraction of blue star forming galaxies in galaxy 
clusters increases with redshifts u p to z ~ 1 (Butcher- 
Oemler effect, iButcher fc Oemlerl Il978l I1984D . Recent 
observations have revealed that there are many star- 
burst galaxies even i n high density reg i ons at z > 1.5 
(lHavashi et all 120101: Iffilton et alj|2010t iPapovich et all 
I2010t iTran et all 120101: iFassbender et all 120111 ). How- 
ever, as their color-magnitude diagrams show a promi- 
nent red sequence, a large fraction of passive galaxies 
also exist in these clusters, and thus formation epoch 
of su ch massive ga l axies seems to be much earlier. In- 
deed, iGobat et al.l (|2011[) have found that there are al- 
ready mature galaxies in a cluster at z = 2.07 which is 
the most distant X-ray cluster to date. It is essential 
to survey dense regions at higher redshifts of z > 2 in 
more detail to investigate the site where a large fraction 
of progenitors of local elliptical are evolving vigorously 
and reveal critical processes for their evolution, although 
these recent studies obviously suggest that star forming 
activity of galaxies become more active with increasing 
redshifts. The importance of surveys in high density re- 
gions at z > 2 is supported by the fact that the activ- 
ities of galaxies and AGNs have peaks at z =1-3 (i.e., 
Madau et ail ll996: Hop kins fc Beacomll200l lUeda et all 



2003). Moreover, it should be noted that massive galax- 



ies on t he red sequence disa ppear in proto-clusters at 
z ~ 3 (jKodama et al.l 120071 ) . while such red galaxies 
are already in place in proto-clusters at z ~ 2 (see also 
i Kaiisawa et all l2006al : iKriek et al.l [2008; Do hertv et all 
I2OIOD . 

Proto-clusters at z > 2 are ideal targets to investi- 
gate the environmental dependence of galaxy proper- 
ties at high redshift. High-z radio galaxies (HzRGs) 
are used as good landmarks to search for proto-clusters, 
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because they are thought to be progenitors of massive 
elliptical gala xies located in the center of local galaxy 
clusters (e.g.. iMcLure et aTl 11999ft . Indeed, some over- 
density regions around HzRGs of various galaxy pop- 
ulations such as Lyman a emitters (LAEs), Lyman 
break galaxies (LBGs), distant re d galaxies (DRGs), are 
identified as prot o -clusters (e.g.. iPentericci et al.l 119971: 
Milev et all 12004 Kurk et all l2004bllal; IKaiisawa et al.l 



2006al IKodama et all 120071: iKuiper et al 1 120101 ) . There- 
fore narrow-band imaging surveys targeting Ha emitters 
(HAEs) around the HzRG is also an effective method 
to investigate star formation activity in a high-z proto- 
cluster. 

There are already several surveys of Ha emission lines 
in proto-cluster s at z > 2; PKS 1 3 38-262 proto - cluste r 
at z = 2.16 (|Kurk et all I2004U lHatch et al.l 1 201 lh . 
4C+10.48 proto-cluster at z = 2.35 (lHatch et al.1 2011 ). 
4C23.56 pro t o-clus ter at z = 2.48 (jTanaka et alJboill i. 
lHatch et alj (|2011l ) have found statistical excesses in 
number density of Ha emitting galaxies in the vicinity 
of the radio galaxies of PKS 1338-262 and 4C+10.48 by 
a factor of > 10 compared with that of general fields. The 
Ha luminosity functions in the proto-clusters have simi- 
lar shapes to those in the general fields, but the normal- 
ization of the luminosity function in the pr oto-clusters is 
a facto r of 13 higher than that in the fields. Tana ka et al.1 
<|201 If ) have found a statistical excess in the number den- 
sity of Ha emitters in the 4C23.56 proto-cluster by a 
factor of 5 compared to a general field. Combining with 
mid-infrared photometric data which traces reradiation 
by dust, they revealed that active star formation must 
been occurring in the proto-clusters, which is as active 
as that in the general fields at similar redshifts. They also 
suggest that it is probable that star formation activity 
in proto-clusters gets stronger towards higher redshifts. 
Therefore, it is interesting to probe star formation activ- 
ity in higher redshift proto-clusters to confirm this trend. 

In this paper, we present results of our Ha emission 
line survey in the 4'x7' region around a radio galaxy 
USS 1558-003 at z = 2.53 with the NB2315 narrow-band 
filter (A c = 2.313/xm,AA = 0.027Aiin). This filter is de- 
signed to perfectly match to this particular proto-cluster 
(Figure [T]). The redshift of 2.53 is the highest where 
Ha emission lines can be traced by near-infrared imag- 
ing on a ground-based tel escope. IKaiisawa et al.l (|2006a) 
and IKodama et al. (2003) have reported a statistical ex- 
cess of bright DRGs around this radio galaxy as an evi- 
dence for the existence of a proto-cluster associated to 
the radio galaxy at z = 2.53. In the current paper, 
our studies have been conducted under the MAHALO- 
Subaru project (MApping H- Alpha and Lines of Oxygen 
with Subaru; Kodama et al., in prep). This project aims 
to reveal the environmental dependence of star forming 
activities by wide-field narrow-band imaging of Ha or 
[O n] emission lines in (proto-) clusters and general fields 
at 1.5 < z < 2.5. Refer also to Kodama et al. (in prep) 
for the details, such as the method of selection of emis- 
sion line galaxies and the estimation of galaxy proper- 
ties (star formation rates and stellar masses). Our pre- 
vious studies have demonstrated that the narrow-band 
imaging is very powerful and useful to completely sam- 
ple star forming gal axies with emission l i nes down to a 
certain lim i ting flux (IKodama et al.ll2004l: iKovama et al.1 
IM0ll23lll[Havashi et al.ll20 1 011201 lb iTadaki et al.ll2011t 
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Fig. 1. — (a) The broken line shows a response curve of MOIRCS 
K s broad-band filter, and the solid line shows that of MOIRCS 
NB2315 narrow-band filter (A c = 2.313^m,AA = 0.027Atm). (b) 
The response curve of NB2315 is shown as a function of redshift. 
The relative line-of-sight velocity with respect to the radio galaxy 
redshift (z = 2.53) is also shown for the case of Ha emission line. 



iTanaka et al~ll2011l ). 

The structure of this paper is the following. Obser- 
vations and available data are described in § Our 
sample of Ha emitters at z — 2.53 around the proto- 
cluster are selected from the photometric catalog in § [3j 
We investigate the spatial distribution and colors of the 
Ha emitters in § |4] Discussions on star forming ac- 
tivity, properties of Ha emitters on the red sequence, 
and the Ha emission of the radio galaxy will follow in 
§ [SJ Finally, we summarize our results of this paper in 
§ [6l Throughout this paper, magnitudes are presented 
in the AB system, and we adopt cosmological parame- 
ters of h = 0.7, fi TO = 0.3 and = 0.7. Vega mag- 
nitudes in J, H and K s , if preferred, can be obtained 
from our AB magnitudes using the following calibra- 
tions: J(Vega)=J(AB)-0.94, (Vega)=iJ(AB)-1.38, 
and iv"(Vega)=iv"(AB)-1.86, respectively. At z = 2.53, 
1 arcmin corresponds to 0.483 Mpc (physical) and 1.705 
Mpc (comoving), respectively. 

2. OBSERVATION AND DATA 

Observations of the USS 1558-003 proto-cluster were 
conducted on 2011 March 11, April 17 and 29 as a Sub- 
aru open-use intensive program (S10B-028, PI: T. Ko- 
dama). The optical and near- infrared (NIR) images 
were taken with Subaru Pr ime Focus Camera (Suprime- 
Cam; iMivazaki etaLl 120021) and Multi-O bject Infra-Red 
Cam e ra and Spectrogra ph (MOIRCS; llchikawa et al.l 
l2006t I Suzuki et al.l 120081) on Subaru telescope, respec- 
tively. Suprime-Cam has a ficld-of-view (FoV) of 27x34 
arcmin 2 , while MOIRCS has a FoV of 4x7 arcmin 2 . The 
K s and NB2315 data were obtained with two pointings 
with an offset of 1' in right ascension and 1' in declina- 
tion, to neatly cover the dense clumps of Ha emitters 
(see fc|3.2[) which were recognized during the course of ob- 
serving runs. We call these two pointings as Fl and F2. 
Although the panoramic Suprime-Cam image covers the 
entire region of both Fl and F2, the MOIRCS J and H 
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TABLE 1 

Summary of the optical and near-infrared images. 



Filter 


FoV a 


Integration time 


Limiting mag. b 


Seeing 


Instrument 


Observation date 






(minutes) 


(5a) 


(arcsec) 







H 



F1+F2 
F1+F2 
F1+F2 

F2 

F2 
F1+F2 



NB2315 F1+F2 



(Fl 



80 
90 
55 
75 
45 
57 
32, F2: 
203 



25) (Fl: 



27.16 
26.87 
25.75 
24.18 
23.51 
23.65 
23.46, F2: 
23.01 



23.17) (Fl 



0.70 
0.63 
0.66 
0.42 
0.47 
0.66 
0.66, F2: 
0.53 



0.40) 



Suprime-Cam 
Suprime-Cam 
Suprime-Cam 
MOIRCS 
MOIRCS 
MOIRCS 



2011 April 29 

2011 April 29 

2011 April 29 
2011 March 11 
2011 March 11 
2011 March 11, 2011 April 17 



(Fl: 133, F2: 70) (Fl: 22.74, F2: 22.35) (Fl: 0.53, F2: 0.36) 



MOIRCS 2011 March 11, 2011 April 17 



Note. — Finally, the FWHMs of PSF in all the images arc matched to 0.66", except for the S-band image which has a FWHM of 0.70". 
a The pointings of Fl and F2 for K s and NB2315 images have an offset of 1' to the west and 1' to the south, 
k The limiting magnitudes are measured with a 1.5" diameter aperture. 



data are available only in the F2 pointing. Therefore, it 
should be noted that the area used in this study is lim- 
ited to the F2 region of about 4x7 arcmin 2 where both 
optical and NIR data are available. The weather was 
fine during the observations, and the sky condition was 
photometric. The seeing was less than 0.66" in all the 
images except for 0.70" in the B-band. The total integra- 
tion times in the optical and NIR broad-bands range from 
45 to 90 minutes. We took significantly longer integra- 
tion (203 minutes) at the narrow-band NB2315. Conse- 
quently, the available data set consists of six broad-band 
data, B,r',z' in optical and J,H,K S in NIR, and the 
narrow-band data with the NB2315 filter. 

All of the optical and NIR data are reduced in a 
standard manner using a data reduction packag e for 
Suprime-Cam (SDFRED vcr.2.0: iQuchi et al.ll2004l) and 
MOI RCS (MCSREER by I . Tana ka), respectively. Re- 
fer to lHavashi et all po'lOL 1201 1[ ) for the details of the 
reduction procedures. PSFs in all the reduced images, 
except for the _B-band image, are matched to 0.66 arcsec. 
The zero-points of magnitudes are determined using the 
standard stars, GD153 and LDS749B for optical data, 
and FS27, FS137 and G191-B2B for NIR data. The 5a 
limiting magnitudes are 27.16, 26.87, 25.75, 24.18, 23.51, 
23.65, and 23.01 in B, r', z', J, H, K s , and NB2315, 
respectively, which are measured with a 1.5" diameter 
aperture. The details of the observation and data are 
summarized in Table [1] 

3. Ha EMITTERS IN THE USS 1558-003 PROTO-CLUSTER 

3.1. Photometric catalog 

Source detection is performed on the original NB2315 
image with a FWHM of PSF of 0.53"rather than the 
PSF matched image, u sing SExtractor (ver. 2.5.0: 
iBertin fc Arnoutsl [l996). and photometry on all the im- 
ages is conducted by the double-image mode of SEX- 
TRACTOR. Color indices are defined with aperture mag- 
nitudes MAG_APER, which are measured with a 1.5"- 
diamcter aperture, while the total magnitudes are de- 
fined with MAG_AUTO magnitudes which are m easured 
with an elliptical aperture (see also iKronl 119801) . Mag- 
nitude errors are estimated from lcr sky noise at each 
object position taking account of slightly different expo- 
sure times and sensitivities, which is thus corrected for 
the difference in depth over the entire FoV. The regions 

5 http://www.naoj .org/staf f /ichi/MCSRED/mc sred.html 



near the edges of each image where the exposure time 
is shorter than 75% of the total in each combined image 
are masked. If an object is not detected in a broad- 
band filter, a 2a limiting magnitude is used to estimate 
the upper /lower limits in color indices or the upper lim- 
its in magnitudes. We make sure that the stellar colors 
are consistent w i th th ose of stellar templates given in 
iGunn fc Strykerl |l983) in order to check the zero-points 
of magnitudes. From such comparison, the zero-point 
magnitudes in B, r' , z' and J are corrected by ~ 0.15 
magnitude at most so that the stellar colors are in good 
agreement with those of the stellar atlas. Moreover, 
magnitudes are corrected for the Galactic absorption 
by the following magnitudes; A(B)=0.62, A(r')=0.41, 
A(O=0.23, A(J)=0.14, A(-ff)=0.09, A(iir 8 )=0.06, and 
A(NB 912)=0.05 which are derived from the extinction 
law of iCardelli et al.l ((1989) on an assumption o f Ry = 
3.1 and E(B-V) = 0.155 based on lSchlegel eTaTl (fl998h . 

As a result, 1,035 objects are detected in the NB2315 
image at more than 5a level (i.e., 23.01 mag. in the 
F1+F2 region). Among them, 754 galaxies are distin- 
guished from 281 stars based on their B — z' and z' — K s 
colors. This method to s e parat e galaxies from stars 
are devised by iDaddi et al.l (|200l ( see also iKong et al.1 
l200l . 

The detection completeness is investigated as follow^. 
First, artificial objects with given magnitudes and Gaus- 
sian profiles with FWHM of 0.53" are distributed on the 
NB2315 image. Source detection is conducted in the 
same manner as described above, and then the fraction 
of detected artificial objects is calculated. The detection 
completeness is found to be more than 74 (83) % down 
to 23.0 (22.0) magnitude in NB2315. 

3.2. Selection of Ha emitters 

Galaxies with a nebular emission line that happens to 
enter the NB2315 narrow-band filter should be observed 
much brighter in NB2315 magnitude compared with the 
K s broad-band magnitude which samples primarily the 
underlying continuum flux density. We thus apply the 
following criteria to select such galaxies with a emission 

6 This completeness refers to point sources, and the one for 
extended objects can be relatively shallower than that. However, 
the sizes of galaxies at z ~ 2.5 are not generally large enough to be 
resolved in the current data, and the assumption of point sources 
may not be so impractical. 
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Fig. 2. — The color-magnitude diagram of K s — NB2315 ver- 
sus NB2315. The solid curves show the boundary of 3<r excess 
given by the equation {TJ. The broken line shows the 2cr limit in 
K s — NB2315 color. The black dots located above the solid curve 
indicate the NB2315 emitters. Some objects with large enough 
color excesses are not classified as emitters due to the shallow data 
at the individual object positions. Blue open circles with the black 
dots represent the Ha emitters associated to the radio galaxy at 
z = 2.53, which are identified by the color selection criteria in r' — J 
and J — K s (equation (|3j and Figure [3} . 

line at ~ 2.313/mi; 

K s - NB > -2.51og(l -y/fl mK3 + fL NB /f NB ) + 0.1, 

(1) 

K S -NB> 0.35, (2) 

where f$ a is the 3cr sky noise flux in each band and /atb 
is the NB2315-band flux density (Figure EJ. The con- 
tinuum flux density is corrected for a color term by O.f 
magnitude as given in the right most side of equation 
(Q}, because the effective wavelength of K s filter is dif- 
ferent by 0.163/zm from that of NB23f 5 filter (Figure CD) 
We no te that the i^ s -NB2315 colors of I Coleman et all 
(1980) templates redshifted to z = 2.53 are actually dis- 
tributed around 0.1 magnitude. The second criterion 
of X S -NB>0.35 corresponds to the observed equivalent 
width larger than 79A (22A in the rest-frame if a galaxy 
is located at z=2.53) after being corrected for the color 
term of 0.1 magnitude. We select galaxies with color 
excesses in if s -NB2315 greater than 3tr photometric er- 
rors. This means that a emission line with a flux larger 
than 2.6xl0 -17 erg s" 1 cm~ 2 can be firmly detected. 
If the color excess is due to a Ha+fNn] emission line 
pair at z = 2.53, the limiting line flux corresponds to 
L(Ha)=l. l x 10 42 erg s' 1 and a dust-free SFR of 8.6 
Mq yr" 1 (|Kennicuttlll998[ ). where t he ratio of [Nil] t o 
Ha([Nn]/Ha) is assumed to be 0.22 (|Sobral et alll2012D . 
We also find 12 galaxies whose line fluxes are larger than 
this limit, but color excesses are smaller than the 3<r level 
due to their location near the edge of the images where 
the depth is shallower. These galaxies are not included 
in our emitter sample. As a result, we select 102 NB2315 
emitters out of 754 galaxies over the ~ 4 x 7 arcmin 2 area 
(Figured]). 

Although we aim to select Ha emitters at z = 2.53, 
it is possible that some of the emission lines detected 
with the NB2315 filter are some other lines at differ- 
ent redshifts such as H/3 or [Om] at z ~ 3.6, [On] at 
z = 5.2, and Paschen series lines at z < 1.0. How- 
ever, it is unlikely that the bulk of the NB2315 emitters 
originate from H/3, [Om], or [On] lines at z > 3.5, be- 
cause the color excesses of the emitters are quite large 



as shown in Figure [5J Note that the large color excesses 
we see in if s -NB2315 may correspond to unrealistically 
large emission line fluxes, if they are all located at such 
high redshifts. Indeed, all but seven NB2315 emitters 
are detected in i?-band and 64 out of 68 Ha emitters 
classified below also have the detection. If a nebular line 
detected by the NB2315 filter would be either an [Om] 
or an [On] line, the line luminosity of [Om] (or [On]) 
is larger than 2.0(4.7) x 10 42 erg s _1 , re spectively. The 
luminosity functions o f [Qui] at 0=0.84 (iLv et all 120071) 
and [On] at z = 1.47 (jSobral et al.ll2012D imply that no 
more than one [O m] or [O n] emitter with such a high lu- 
minosity can be included in our NB2315 emitter sample 
with limited field coverage, if any. Therefore, we believe 
that there is few contamination from the lines at z > 3.5 
in our NB2315 emitter sample. Moreover, it should be 
noted that our NB2315 survey is targeting the plausible 
proto-cluster at z=2.53 associated to the radio galaxy, 
and we must be preferentially detecting Ha emitters at 
z = 2.53 rather than other high -z emitters in the ge neral 
field. However, as reported by iGeach et all (j2008[ ) who 
conducted the H2SI narrow-band (A c = 2.121/im) imag- 
ing with WFCAM on UKIRT to search for Ha emitters 
at z = 2.23 in the COSMOS field, other lower- z emission 
lines such as Paa and Pa/3 at z < 1 can be significant 
contaminations. In order to discriminate Ha from other 
lower- z contaminant lines, we set the following color se- 
lection criteria on the color-color diagram, r' — J versus 
J-K s ; 

J — K s > 0.7 (r' - J) + 0.05 for r' — J < 1.9 

J -K s > 1.38 for r' - J > 1.9. [6) 

This classification with r'JK s colors is an analogue of 
the BzK c olor selection of ga laxies at 1.4 < z < 2.5 
devised by iDaddi et all (|2004[ ). but it is optimized to 
identify galaxies at z > 2.5. Figure [3] shows the color 
tracks of model spectra with various star formation his- 
tories and fixed formation redshift of Zf = 5 based on 
Bruzual fc Charlot| (|2003l) . and four SED templates in 
iColeman et all (|1980( ) which are redshifted from z = 0.0 
to z — 4.0. This figure suggests that the criteria enable 
us not only to select DRGs as galaxies with red colors 
(J — if s (vega)> 2.3), but also to isolate star forming- 
galaxies at z > 2.5 from those at lower redshifts. Al- 
though no dust extinction is taken into account for the 
color tracks, the arrow in the figure sho ws the redden- 
ing ve ctor of E(B — V) = 0.2 based on the lCalzetti et all 
( 2000) extinction curve. This indicates that our selection 
criteria are relatively free from the dust reddening effect, 
as the vector is almost parallel to the slanted boundary 
line. Indeed, a group of emission line galaxies which are 
likely to be Ha emitters at z — 2.53 seem to be located 
separately from the bulk of galaxies at lower redshifts as 
shown in the right panel of Figure [3] Thus, our selection 
criteria work well to select Ha emitters at z = 2.53. 

We attempted to apply other color sel ection meth- 
ods s uch as the original Bz'K s selectio n (Dadd iet all 
2004f) and JHK f sel ection as well (jKaiisawa et all l2006a; 
Kodama et al.1120071 ). However, the Bz'K s selection is ef- 
fective to select galaxies at 1.4 < z < 2.5, while we aim 
to pick out Ha emitters at z ~ 2.5 among the NB2315 
emitters. Since they are located just at the upper edge 
of the redshift range of the Bz'K selection, the selection 
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Fig. 3. — (Left panel) The color-color diagram of r' — J versus J — Ks. Red, orange, and gre en solid lines show color t racks of model 
SEDs for three types of star formation histories with fixed formation redshift of Zf = 5 based on Bruzual & Chariot (2003). Dotted lines 
show the color tracks created by four empirical spectral templates from elliptical to irreg ular type ga laxies compiled by ICoIeman et alj 
(1980). Dust extinction is not corrected, but the absorption by neutral hydrogen given by Madau (1995) is tak en in t o acco unt. The arrow 
shows a reddening vector o f E(B — V) = 0.2, which is estimated from the dust extinction curve of Calzetti ct al. (2000). The numbers 
shown along the tracks for Bruzual & Chariot (2003) model indicate rcdshifts, and the open squares correspond to z = 2.5. The black 
folded solid line shows the boundaries of our selection criteria to identify Ha emitters among the NB2315 emitters (equation (£3j ) . (Right 
panel) Same as the left panel, but the observed galaxy colors are plotted. Gray dots show the NB2315-detected galaxies. Blue open circles 
show our sample of Ha emitters, and black open squares show the other contaminant line emitters probably located at lower-z (z < 1). 



of Ha emitters may well be quite incomplete. On the 
other hand, it is found that the JHK S selection do not 
work well due to the shallowness of .ff-band data and 
the incompleteness in the selection of blue star-forming 
galaxies. Therefore, we judge that the classification with 
r' JK S colors works best among them to identify Ha emit- 
ters at z = 2.53. Instead of the color selection, one may 
consider the use of photometric redshifts to discriminate 
among different emission lines at different redshifts. We 
also derive photometric redshifts of our NB2315-detected 
galaxies using six ba nds photometry of B r' , z' , J, H, K s 
and the EAZY code (|Brammer et al.f 2008). However, we 
find that it is difficult to cleanly discriminate Ha emit- 
ters from other possible lines, especially those at z > 1.0, 
based on the photometric redshifts. In the end, we decide 
to use the r'JK s diagram to identify Ha emitters. Con- 
sequently, we select 68 Ha emitters in total (Figure [3]). 

3.3. Selection of Distant Red Galaxies 

DRG is a class of galaxy populations at 2 < z < 4 
wit h red colors of J — K s > 1.38 or J — K s (vega,)> 
2.3 (jFranx et al.ll2003T ) which consist of dusty starburst 
galaxies and quiescent galaxies. To select such red galax- 
ies in this region, the DRG color selection is applied to 
the i^g-selected catalog which is made with the same 
procedures as those described in § 13. II except that source 
detection is performed on the K s image. The iC s -selected 
catalog is used only for the selection of DRGs, since such 
catalog enables us to select DRGs more completely down 
to fainter magnitude than the NB2315-limiting magni- 
tude. The catalog includes 1,340 objects brighter than 
5cr limiting magnitude in K s . Among them, 1,002 objects 
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Fig. 4. — Differential number densities of DRGs and ii's-selected 
galaxies around the USS 1558-003 proto-cluster are shown in up- 
per panel, and cumulative number densities of DRGs are shown in 
lower panel. Filled circles show the number densities corrected for 
detection incompleteness, while open circles show those without 
corrections. Filled squares and tr iangles show the resu lts obtained 
by MODS deep an d wide surveys (Kajisawa ct al. 2011), and aster- 
isks show those of Kaiisawa ct al. (2006b). The results of MODS 
deep and wide are not corrected for detection incompleteness. Er- 
ror bars are estimated based on Poisson statistics. Vertical broken 
line shows the 5<r limiting magnitude in K s . 
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TABLE 2 

The number and number density of Ha emitters (HAEs) 
and DRGs. 



Region 



Area 
(arcmin 2 ) 



Number 



Density 
(arcmin -2 ) 







HAE 


DRG 


HAE 


DRG 


Clump 1 


3.36 


15 


12 


4.46±1.15 


3.57±1.03 


Clump 2 


1.64 


20 


8 


12.2±2.73 


4.88±1.72 


Clump 3 


0.94 


8 


3 


8.51±3.01 


3.19±1.84 


All clumps 


5.94 


43 


23 


7.24±1.10 


3.87±0.81 


Others 


21.16 


25 


19 


1.18±0.24 


0.90±0.21 


Entire field 


27.10 


68 


42 


2.51±0.30 


1.55±0.24 



Note. — Errors in the number density arc estimated based on Pois- 
son statistics. 

are classified as galaxies, while the remaining 338 objects 
are classified as stars. Consequently, 42 DRGs are se- 
lected from the if s -selected catalog. It is expected that 
red galaxies associated to the proto-cluster at z = 2.53 
tend to be seen as DRGs, although DRGs in general can 
be located in a much wider redshift range of 2 < z < 4. 

Figure [4] shows the number densities of DRGs and K s - 
selected galaxies which are corrected for detection in- 
completeness estimated in the same manner as in ij3.ll 
For comparison with other studies, we also plot the re- 
sults of the gene ral blank field survey, MOIRCS Deep 
Survey (MODS) (jKaiisawa et al.ll2006bi |20H where the 
deep and wide surveys have FoVs of 28.2 and 103.3 
arcmin 2 , respectively. The number densities of DRGs 
as well as if s -selected galaxies in the USS1558 proto- 
cluster do not show strong overdensities compared to 
those in the MODS, although so me excesses by a factor 
of a few are seen at several bins. iKodama et al.l ([2007D 
found that this region shows an over density of bright 
DRG s compared to the GOODS-S field (Giavalisc o et al.l 
2004) . Figure 0] indeed shows that there is an excess of 
bright DRGs with K s < 21.5 in the USS 1558-003 proto- 
cluster which are almost absent in the general blank field. 
The difference in number densities of DRGs between the 
proto-cluster and the general field is clearer in the cumu- 
lative densities, showing the excess at ~ 2cr level in the 
bin with K s < 21.5. Moreover, the spatial distribution 
of DRGs is not homogeneous but clustered around the 
radio galaxy (see N4.1|) . The fact that there is no signifi- 
cant excess of faint DRGs in the proto-cluster region may 
suggest that most of the faint member galaxies in this 
proto-cluster have bluer colors due probably to on-going 
star formation activities, and that passive red galaxies 
would be gradually emerging later on. Although we do 
not have yet any spectroscopic confirmation of member- 
ship of the DRGs except for the radio galaxy itself, it 
is likely that the region around the USS 1558-003 radio 
galaxy is a proto-cluster at z — 2.53. 



4. RESULTS 



4.1. 



Spatial distribution of Ha emitters 

Figure [5] shows the spatial distribution of our 68 Ha 
emitter (HAE) candidates associated to the USS 1558- 
003 radio galaxy at z = 2.53 located at (0,0) in the figure. 
The 42 DRGs are also plotted with red filled circles. We 
find that there are three outstanding regions where HAEs 
and/or DRGs are strongly clustered, which we hereafter 
call clump- 1, clump-2, and clump-3. In fact, as sum- 
marized in Table [2j the number densities of HAEs and 



DRGs in these three clumps are all higher than those of 
averaged values across the observed fields by factors of 
2-5 for HAEs and 2-3 for DRGs. The clump-1 is the 
vicinity of the radio galaxy, and it contains both HAEs 
and DRGs around the radio galaxy. If we assume that 
the central dominant radio galaxy grows to a cD galaxy 
in the future, this region may correspond to the cen- 
tral part of the cluster. The fact that DRGs are also 
clustered in this clump may also suggest that this part 
of the proto-cluster is the oldest. The clump-2 is the 
most conspicuous, densest association of the HAEs and 
DRGs, and is located at ~ 3.2' (about 1.5 Mpc in physi- 
cal scale) away from the radio galaxy to the south-west. 
The clump-3 is a smaller group of Ha emitters located 
in between the clump-1 and the clump-2. 

It is interesting to note that the Ha emitters are 
more strongly clustered towards the south-west clumps 
(clump-2 and clump-3) rather than in the immediate 
surrounding region around the radio galaxy (clump-1). 
They constitute a part of large scale structure hosting 
the radio galaxy and clump-1, and they would all merge 
together in the near future to form a more massive single 
cluster around the radio galaxy. It is obvious that this 
proto-cluster region is not relaxed yet and just in the 
process of galaxy assembly from the the surrounding re- 
gions. We note that such distribution of t he Ha emitters 
is sim ilar to that of DRGs reported by Kodama et al. 
;2(I(I7) based on NTT/SOFI imaging data. In fact, we 
also confirm that the DRGs tend to be located in and 
along the structures traced by the Ha emitters based on 
our deeper MOIRCS data (Figure [5]). Furthermore, we 
find that these DRGs also meet the BzK criteria which 
select galaxies primarily at 1.4 < z < 2.5, suggesting that 
significant fraction of these DRGs are likely to be phys- 
ically associated to the proto-cluster hosting the radio 
galaxy at z = 2.53. 

The three clumps in the 1558-003 proto-cluster host a 
large number of Ha emitters. The surface number den- 
sity of Ha emitters in the three clumps is ~37±13 times 
larger than those in the redshift slice at z = 2.2 in the 
GOODS-North fiel d (Ha emission sur vey with NB209 
narrow-band filter; [Tadaki et all 1201 ID . Similarly, we 
also find that it seems that the surface number density is 
~14±3 and 17±4 times larger than those in the redshift 
slices of z =2.2 and 2.53 in the SXDS field (Ha emission 
s urvey with NB209 ar id NB2315; Tadaki et al. in prep.). 

iTanaka et al.l (|2011f ) conducted a Ha emitters survey 
in the field around the 4C 23.56 radio galaxy at z = 2.48 
with CO narrow-band filter on MOIRCS, and found 11 
Ha emitters candidates to flux down to ~ 7.5 x 10~ 17 
erg s _1 cm -1 and rest-frame EW > 50A over a 23.6 
arcmin 2 area, which is similar to our survey area. It is 
found that Ha emitters are distributed on the east side 
of the radio galaxy 4C 23.56, and there is a clump of 
Ha emitters ~ 2 Mpc (in comoving scale) away from the 
radio galaxy. Such an offset distribution of Ha emitters 
from the radio galaxy is similar to that in our USS 1558- 
003 field, although our Ha emitters are more strongly 
clustered than the 4C 23.56 field. However, note that 
the EW cut and the limiting flux used for Ha emitter 
selection are s lightl y different between our survey and 
ITanaka et all (|20ll . Our survey enables us to sample 
Ha emitters with fainter line fluxes and smaller EWs. 
If we apply the same EW cut and the limiting flux as 
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in iTanaka et all (|2011[ ) to our sample in the USS 1558- 
003 field, the number of Ha emitters reduces to 27, but 
it is still considerably larger than that in the 4C 23.56 
field. The discovery of clumps of Ha emitters in these 
proto-clusters around the radio galaxy at z ~ 2.5 clearly 
indicates that we are witnessing the process of mass as- 
sembly of clusters at their early stage when galaxies are 
vigorously assembling to form dense cluster cores while 
they are acti vely forming stars. 



Moreover, lHatch et al 



(|2011l ) report the studies for 
two proto-clusters around the radio galaxies MRC 1138- 
262 at z = 2.16 and 4C+10.48 at z = 2.35. These stud- 
ies suggest that proto-clusters at z > 2 tend to show 
the statistical excess of Ha emitters compared to the 
general fields at similar redshifts. However, they found 
that Ha emitters in the 4C+10.48 region are not strongly 
clustered, wh ile those in the M RC 1138-262 region are 
clustered fsee lKurk et alJl2004b[ ) . Moreover, both proto- 
clusters show inhomogeneous spatial distribution of Ha 
emitters around each radio galaxy, and most of the Ha 
emitters are distributed only on one side of the radio 
galaxies, similar to 4C 23.56 (jTanaka et al.l 1 2 llh and 
USS 1558-003 (this study). 

Another interesting result we find is that the radio 
galaxy show an extremely extended Ha emission spa- 
tially as shown in Figure[BJ In the NB2315 image, we can 
notice that the structure of the Ha emission is stretched 
in the northeast-southwest direction. The size of the Ha 
emission is ~ 4.5"which corresponds to ^36 kpc in phys- 
ical scale. We discuss the extended Ha emission of radio 
galaxy itself in i j5.2l 

4.2. Color-magnitude diagram 

Color-magnitude diagram is a powerful tool to inves- 
tigate the properties of galaxies in clusters or proto- 
clusters. It is well-known that galaxy clusters at low red- 
shifts are dominated by quiescent galaxies which make 
up the red sequence on the color-magnitude diagram. 
The tight sequence of red galaxies is one of the promi- 
nent features seen in galaxy clusters. However, when 
and how do such red quiescent galaxies form and then 
evolve? The answer to this question still remains un- 
clear, although there are some evidences suggesting that 
the blue star forming galaxies in the early Universe 
evolve and become red quiescent galaxies in high den- 
sity regions during the redshift interval o f 2-3 (e.g., 
Kaiisawa et al.ll2006al IKodama et"aT1l2007t | Kriek et al.l 
20081: IDohertv et al.l 120101 : IGobat et alTEnl . To ad- 
dress this, it is essential to investigate further the color- 
magnitude diagram for galaxy cluster and proto-cluster 
at high redshifts when the clusters are vigorously evolv- 
ing. 

Figure [7] shows a color-magnitude diagram of J — K s 
versus K s for the galaxies in the observed field towards 
the proto-cluster USS 1558-003. The Ha emitters asso- 
ciated to the pr oto-cluster are plotted by open circles. 
As presented by IKodama et al.l (|2007fl , we may be able 
to recognize a sequence of DRGs at J — K s ~ 1.5 which 
corresponds to the color of qui escent galaxies with for- 
mation redshift of 3 < Zf < 4 (jKodama et al.lfl~998h . as 
indicated by the tickmarks in the figure. It is not surpris- 
ing that the majority of the Ha emitters are located on 
the bluer side of the diagram. However, it is intriguing 
that some Ha emitters have very red colors satisfying 
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Fig. 5. — The spatial distribution of Ha emitters. Blue open cir- 
cles show Ha emitters at z = 2.53. Red filled circles shows DRGs, 
and black dots are NB2315-detected galaxies. North is up, and 
east is to the left. The origin of the coordinates is the position of 
the USS 1558-003 radio galaxy. Three gray circles are the regions 
where Ha emitters and DRGs are strongly clustered, which are de- 
fined as clump-1, clump-2 and clump-3, respectively. Two regions 
enclosed by gray solid lines show our MOIRCS pointings (Fl and 
F2). 




Fig. 6. — The close-up views of the USS 1558-003 radio galaxy 
at z = 2.53 in the NB2315, K s and NB2315-,f£'. J images. Note 
that the FWHM values of PSF in these Ks and NB2315 images 
are 0.40" and 0.36", respectively, since these images are created by 
combining only the frames taken under excellent seeing conditions 
~ 0.4". The length of the bar in the NB2315-if s image indicates 
the angular size of 4.5", which corresponds to 36.2 kpc in physi- 
cal scale. The NB2315-ft's image clearly indicates that the radio 
galaxy has an extremely extended Ha emission. 



the DRG criterion and constitute mainly the faint end 
of the red sequence. We also note that there is a sig- 
nificant color scatter among the red Ha emitters. Sim- 
ilarly red Ha emitters and dusty star forming galaxies 
have also been recogni zed by recent observations of clus- 
ters at lowe r redshifts (|Geach et al.l l2006: Ko vama et al.l 
120101 1201 If ). They tend to be dusty star forming galax- 
ies and are considered to be in the transition phase from 
active star forming galaxies to passive quiescent galax- 
ies probably under the influence of some environmen- 
tal effects, because these are preferentially found in the 
medium density regions of galaxy clusters at intermedi- 
ate redshifts (z < 0.8) when and where we see a sharp 
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Fig. 7. — Color-magnitude diagram of J — K a versus Ks. Blue 
open circles show Ha emitters at z = 2.53 and black dots show 
all the galaxies in the observed field. Long-dashed, short-dashed, 
and dotted lines show 5cr, 3cr, and 2<r limits in color and mag- 
nitude. Red solid horizontal line shows the color corresponding 
to J — _R" S =1.38, which is a criterion for selecting DRGs. The 
J — K s colors of the red quiescent galaxies with formation rcdshifts 
of 2/=3, 4, and 5 are shown by the tickmarks at the left edge, which 
are estimated by thejKodama ct al. (1998) model. Dotted curves 
show the iso-stellar mass curves for lxlO 1:L M0and lxl0 10 M^ 
spectively. 
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tran sition in the distribution of galaxy colors from blue to 
red (jTanaka et al.ll200l iKwama et al.|[2008h . Therefore 
it is suggestive that the faint end of the red sequence is 
just being built by those transitional galaxies recognized 
as the red Ha emitters. 

For a galaxy at z = 2.53, its K s band luminosity is still 
a good proxy for the stellar mass. However, the mass-to- 
luminosity ratio (M+l Lk b ) is dependent on SED, and it 
should be corrected to get more precise stellar mass from 
the .fG-band luminosity. We use a simple method to es- 
timate stellar masses of the Ha emitters based on their 
-Kg-band magnitudes and J — K s colors. The relationship 
between J — K s colors and M+/Lk s are approximately 
estimated usin g the s t ellar p opulation synthesis model of 
IKodama et al.l (|1998l 119991) where we constructed a se- 
quence of model with varying bulge-to-disk ratios. We 
note here that because of the SED degeneracy among 
age, metallicity, and dust extinction, the relationship be- 
tween J — K s and M+/ Lk s is relatively insensitive to the 
detailed modeling of stellar populations excep t for the 
effect of IMF variation. We here assume the ISalpeterl 
(fT955l) IMF, and the stellar masses in the models are 
scaled accordingly. From this experiment, we establish 
the following relationship between stellar mass (M*), K s 
magnitude, and J — K s color; 

log 10 (A4/10 n Af o ) = -0A(K s -K u )+Alog w M, (4) 

where K\i is the ii" s -band total magnitude corresponding 
to the stellar mass of 10 11 M Q for a passively evolving 
galaxy formed at z=5 and observed at z=2.53, which 
is estimated to 22.41. The mass-to-luminosity ratio is 
corrected for depending on the J — K a color by, 

Alog 10 M = 0.12 - 1.84 • exp(-1.35 • (J - K s )). (5) 



The iso-stellar mass curves of 1 x 10 11 M Q and 1 x 10 10 M© 
are shown in Figure [71 All the Ha emitters but the radio 
galaxy have stellar masses of < 1O 11 M0. Among them, 
the bluest Ha emitters with J — K s < 0.7 are less massive 
galaxies with stellar masses of < 10 10 Mq. This gives a 
dichotomy in the distribution of star forming galaxies on 
the stellar mass-color plane, separated at ~ 10 10 M Q , be- 
low which the emitters are the bluest hence the youngest 
or the most actively star forming. 

It should be also noted that the brightest Ha emit- 
ter is the radio galaxy, and it is located at the bright- 
est end of the red sequence. Note that similar re- 
sults axe_als^_Jojxnd_Jn other proto-cluster s at z =2 
2.5 (jTanaka et al.l 120111: lHatch et al.l 120111) . Its stel- 
lar mass is estimated to ~ 10 12 M^, which i s typ- 
ical for HzR Gs (e . g.. iRocca-Volmerange et al.l l2004t 
iSevmour et"al]|2007HHatch. et al.ll2009l ). Since the radio 



galaxy hosts an AGN, it is not straightforward to quan- 
tify its star formation rate from the emission line strength 
as it is affected by photoionization by AGN f ij5.2p . 

4.3. Star formation activity 

The Ha luminosity is a good indicator of star formation 
rate in the individual galaxy, because Ha emission line is 
originated from the photoionization by young O and B 
stars, and thus it probes the instantaneous star formation 
activity on a time scale of < 20 Myr. The advantage of 
using Ha line to estimate star formation rate is that it is 
well-calibrated and widely used in the studies of galaxies 
both in the local and high-z Universe. Moreover it is 
relatively insensitive to dust extinction compared to [O n] 
and UV luminosities. 

The Ha luminosity is estimated from the emission line 
flux entering into the narrow-band NB2315, which is 
the combination of Ha and [Nil] lines in most cases. 
Since the emission lines do not enter the broad-band fil- 
ter for the Ha emitters at z = 2.53, as shown in Fig- 
ure [I] the flux densities in the narrow-band and the 
broad-band filters are expressed as Jnb — /continuum + 
F(Ka + [Nn])/AiVB and f BB = /continuum, and thus the 
emission line flux is: 



F(Ha + [Nu]) = {f NB - f BB )-A 



NB, (6) 

where A NB = 271 A is FWHMs of the filters, and 
/continuum is the flux density of continuum level of a spec- 
trum. It should be noted that the fluxes are derived un- 
der the assumption that the emission lines are shifted 
into the center of the narrow-band filter. Since the re- 
sponse curve of the filter is not perfectly top-hat (Fig- 
ure [lj , the fluxes can be underestimated if the emission 
lines are not located in the center of the filter. In that 
case, the SFR derived from the flux can be a lower limit. 
The contribution of [Nil] is removed from the line flux, 
F(Ha+[Nn]) by assuming the relation between the ra- 
tio of [Nil] /Ha and the rest-frame equivalen t width of 
EW (Ha+[Nn]) given bv lSobral et all (12011) . 

We derive SFR from Ha flux using the IKennicuttl 
(1998) relation, where Ha flux is corrected for dust ex- 
tinction using the SFR-dependent calibration given by 
iGarn et"al1 (|2010l ) . The dust extinction correction is im- 
portant even for Ha line which is located in the rest- 
frame optical wavelength (A6563). It is often assumed 
that the amount of dust extinction, A(Ha), is unity in 
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Fig. 8. — (Left panel) Dust extinction un-corrected SFRs of the 68 Ho emitters as a function of stellar mass. Those in clumps-1, 2, and 
3 are marked by blue open circles, cyan open squares, and green open triangles, respectively. (Right panel) Same as the left panel, but the 
SFRs are corrected for dust extinction (se e text for details). The solid line shows a main sequence (MS) of star forming galaxies at z ~ 2 
(Daddi ct al. 2007; Rodighicro ct al. 2011), and the broken line corresponds to the SFRs elevated by a factor of 4 compared to the MS for 
a given stellar mass. 

Mq is the radio galaxy USS 1558-003 and it apparently 
shows a considerably high SFR. However, this value must 
be wrong because its Ha line intensity is severely contam- 
inated by the AGN component. In the right panel of Fig- 
ure |51 we plot the main sequence of st ar forming galax- 
ies at z ~ 2 (|Daddi et al.ll2007t) . In iRodighiero et al.l 
(|201 If) . the galaxies whose SFRs are more than 4 times 
larger than those on the main sequence for a given stellar 
are defined as star-burst galaxies. Following this defini- 
tion, a large fraction of galaxies with < 10 10 Mq are 
starburst galaxies, while the massive galaxies are more 
or less located around the main sequence. This is some- 
what expected by the fact that less massive galaxies show 
the bluest colors on the color-magnitude diagram (Fig- 
ure [7]). In contrast, none of massive galaxies shows such 
bluest colors. It should be noted, however, that the red 
Ha emitters with stellar masses larger than 10 10 M Q are 
subject to large amount of dust attenuation, and thus 
they are actually in dusty starburst phase as we will dis- 
cuss later in ij l5.ll Figure|H]shows specific SFRs of the Ha 
emitters plotted as a function of stellar mass. It indicates 
that specific SFR and stellar mass are anti-correlated as 
seen in many other studies. Massive galaxies have lower 
star formation rates for a given stellar mass compared 
to less massive galaxies, indicating that star forming ac- 
tivities have been somewhat truncated and weakened in 
massive galaxies. 

In Figures[5]and|ni Ha emitters in the three clumps are 
marked with open circles, squares, and triangles. It ap- 
pears that the Ha emitters in the clump-2 (i.e., the south- 
west clump ^1.5 Mpc away from the radio galaxy) tend 
to have higher star formation rates than those in other 
regions. To investigate any environmental dependence 
in the properties of Ha emitters more quantitatively, we 
show in Figure [TU] the differential histograms of SFRs, 
specific SFRs, and stellar masses of the Ha emitters di- 
vided by the regions in the left panels, and their cumu- 
lative and normalized number counts in the right panels. 
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Fig. 9. — Specific SFRs of the Ha emitters as a function of stellar 
mass. Dust extinction is corrected (see text for details). Symbols 
are the same as in Figure [8] Three dotted lines correspond to the 
constant SFRs of 10, 100, and 1000 M yr -1 , respectively. 



many papers. However, we adopt the iGarn et al.l ([2010D 
relation to more realistically estimate the intrinsic Ha lu- 
minosity. The specific SFR is also calculated by dividing 
the SFR by the stellar mass estimated in H4.2I 

Figure El shows SFRs of the Ha emitters as a func- 
tion of stellar mass, where dust-uncorrected and dust- 
corrected SFRs (SFR corr ) are shown in separate panels. 
There is a correlation that more massive galaxies tend to 
have higher SFRs, but the dependence is weak because 
there is a significant scatter in SFR for a given stellar 
mass. A large fraction of the Ha emitters have high 
SFR corr larger than 100 M Q yr -1 , meaning that many 
star forming galaxies in this proto-cluster are in the star- 
burst phase. The object with the stellar mass of ~ 10 12 
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Fig. 10. — (Left panels) Distribution of SFR (top), specific SFR (middle) and stellar mass (bottom) for the Ha emitters in the clump-1 
(blue), clump-2 (cyan), clump-3 (green), and the other region (black), respectively. (Right panels) Normalized cumulative number counts 
of the Ha emitters in each region as a function of SFR (top), specific SFR (middle) and stellar mass (bottom). 



There is no large difference in the distribution of all 
these physical quantities among different regions in gen- 
eral, except for the weak tendency that SFRs in the 
clump-2 may be slightly higher. In fact, the median SFR 
of HAE in the clump-2 is 152 M Q yr -1 which compared 
to 97 Mq yr -1 in the clump-1, 99 M Q yr -1 in the clump- 
3 and 106 M© yr -1 in the other regions. However, the 
Kolmogorov-Smirnov (KS) test that is applied on the top 
right panel suggests that the difference does not show 
high enough significance in statistics. The probability to 
accept the hypothesis that the SFRs in these clumps have 
different distributions is less than 2cr at most. Therefore 
we may conclude that star formation activity at this high 
redshift (~2.5) is high everywhere regardless of environ- 
ment, not only in the low density outskirts but also in 
the dense clumps of the proto-cluster. Because we are 
targeting a proto-cluster and its surrounding region, it 
is yet unknown whether such absence of environmental 
dependence at z ~ 2.5 is universal from a proto-cluster 
to another, as well as in the general low density fields 
which are not in the vicinity of proto-clusters. 

Furthermore, it is interesting to know how much star 
formation activity and stellar mass are confined in this 
proto-cluster region as a whole, as most of the structures 
and the Ha emitters are likely to assemble to form a sin- 
gle richer cluster at later times. The integrated stellar 
mass and SFR of the Ha emitters in each clump and 
those in the entire region are shown in Table El The 
uncertainty in the integrated stellar mass and SFR orig- 
inated from photometric errors are small enough not to 
give a large influence on the discussions below. Also, in 
this analysis, we exclude the most massive radio galaxy 
in the clump-1 (~lxl0 12 M ), since its Ha line flux 
is severely contaminated by the AGN component. The 
integrated stellar mass of DRGs associated to the proto- 
cluster is derived by subtracting the component of DRGs 
outside the cluster but projected along the line of slight. 
Here we use the nu mber count of DRGs in the M ODS 
general field survey (|Kaiisawa et al.ll2006bl lUIl . We 
note that the total stellar masses estimated here should 



be taken as the lower limits, because we miss less massive 
galaxies fainter than our observational limits. All clumps 
but clump-3 have a total stellar mass of ~ 10 12 M Q , while 
clump-3 has a factor of ~10 lower mass of ~ lO 11 ^!©. 
The integrated specific SFR of the clump-1 is slightly 
lower than those of the clump-2 and the clump-3, by a 
factor of 1.6 and 3.9, respectively, which might suggest 
that the star formation activity per unit stellar mass is 
slightly lower in the closest vicinity of the radio galaxy 
which may grow to the center of the cluster in the future, 
and the age of the system is slightly older than the other 
clumps. However, the statistics is again too poor for us 
to say anything conclusive yet. 

5. DISCUSSIONS 

5.1. red Ha emitters 

As described in § 14. 2\ there are some H AEs with colors 
redder than J — A' s =1.38. Figure [Til shows the distri- 
bution of the HAEs and DRGs, where the red HAEs 
with J — K s > 1.38 are marked by red open circles. It 
seems that the red HAEs are clustered around the radio 
galaxy (clump-1) and in the south-west clumps (clump-2 
and clump-3). Such concentration of the red emitters in 
high density regions is qualitatively different from what 
is seen in the lower redshift cluster at z ~ 0.8 where red 
HAEs prefer to medium density regions, as mentioned 
in § 14.21 This may indicate that, in this proto-cluster 
at z ~ 2.5, the transition of galaxies is occurring in the 
densest environment rather than in the outskirts. We 
note that most of the red HAEs are less massive than 
10 n M Q as shown in Figure If the dusty starburst 
phase does not last so long, these red HAEs would stay 
on relatively faint end of the red sequence even after they 
quench star formation. If the DRGs without Ha emis- 
sion lines with stellar masses of a few x lO n M0are really 
proto-cluster members, the active phase of these mas- 
sive quiescent galaxies should have been seen at z > 3. 
Some of their progenitors would be like the popula- 
tions of LBGs or LAEs which are in fact often identi- 
fied in over-dense regions (proto-clusters) at z > 3 (e.g., 
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TABLE 3 

The integrated stellar masses and SFRs of HAEs and DRGs. 



Region 


Stellar 


mass (E Ah) [M Q ] 


SFR (S SFR) [M yr" 1 ] 


Specific SFR (S SFR / £ M*) [yr" 1 ] 




HAE(blue) 


HAE(red) 


DRG 


HAE 


HAE+DRG 


Clump 1 


1.3X10 11 


2.3X10 11 


1.2xl0 12 


1.8xl0 3 


1.4xl0" 9 


Clump 2 


3.9X10 11 


2.7X10 11 


1.2xl0 12 


3.4xl0 3 


2.2xl(T 9 


Clump 3 


8.2xl0 10 


9.1 xlO 10 


l.lxlO 11 


l.OxlO 3 


5.4xl(T 9 


All clumps 


6.1X10 11 


6.0X10 11 


2.2xl0 12 


6.2xl0 3 


2.2xl(T 9 


Others 


4.1X10 11 


1.8X10 11 


9.2X10 11 


3.2xl0 3 


2.4xl0~ 9 


Entire field 


l.OxlO 12 


7.8X10 11 


2.8xl0 12 


9.4xl0 3 


2.4xl0~ 9 



NOTE. — The numbers arc also divided into clumps and to the other region. For HAEs, blue and red ones arc separated 
at J — K s — 1.38 (same as the DRG criterion). For DRGs, the field contamination from foreground/background fields along 
the line of sight is statistically subtracted using the general field data (MODS). In estimation of specific SFR for combined 
samples with HAE and DRG, stellar masses are calculated by summing up those for HAE(bluc) and DRG. The radio galaxy 
in the clump-1 is excluded in this table, because its Ha flux is severely contaminated by the AGN component. 



Steidel et all 119981: IMilev et al.l [20(51 IVenemans et al.l 
2005L 120071 : lOuchT et al. 20051: IKuiper et al.H2011li ~ 

On the other hand, as another possibility, the red 
HAEs can be quiescent galaxies with AGN activity in the 
galaxy core, not dusty starburst galaxies. Even in that 
case, the fact that the red HAEs are clustered in high 
density region is very interesting. This is because the 
considered situation suggests that AGN activity is en- 
hanced in the core regions of the clumps, and thus AGN 
must give a significant influence on the evolution of high- 
z galaxies in high density region. Perhaps, AGN activ- 
ity has something to do with process to quench the star 
formation of cluster galaxies, i.e., so-called AGN feed- 
back. Although it is difficult to reveal which hypotheses 
are true between dust starbursts and quiescent galaxies 
hosting an AGN with only the data now available, in any 
case, there is no doubt that galaxies are activated in the 
densest regions such as clumps at z ~ 2.5. 

In lower-z clusters (z < 1.0), the fraction of star 
forming galaxies to all cluster members decreases as 
we go closer to cluster centers or as the number den- 
sity of galaxies increase s (e.g., iKodama et al.l l2004t 
IKovaina et aTll2010l l20ll . In this USS 1558-003 proto- 
cluster, however, it is difficult to completely sample red 
quiescent galaxies, since DRGs have broad redshift dis- 
tribution such as 2 < z < 4 and also contain both popu- 
lations of quiescent galaxies and dusty starbursts. This 
makes it hard for us to quantify the fraction of star form- 
ing galaxies in this proto-cluster. However, even if we 
assume that all the DRGs are quiescent galaxies asso- 
ciated to the proto-cluster at z = 2.53, a large fraction 
of proto-cluster members are star-formers, and there are 
much fewer red quiescent galaxies. 

5.2. extended Ha emission in radio galaxy 

As shown in Figure El the radio galaxy USS 1558-003 
has quite an extended Ha emission with the apparent 
scale of ~ 4.5" or the physical scale of ^36 kpc. Al- 
though the spatially extended nebular line emission (Lya 
and Ha) ha ve already been recognized around this ra- 
dio g alaxy (|Villar-Martm et all 120071 : iHumphrev et al.l 
2008), we here reveal for the first time the two- 
dimensional distribution of Ha emission. Such extended 
nebular emission is a remarkable featu r e commonly seen 
in H zRGs (e.g., iHeckman et alJH991t IHumphrev et al.l 
2007). In particular, it is well-known that they show ex- 
tended Lya emission aligned with radio structure (align- 



ment effect). Our extended Ha emission is tilted with 
the position angle (PA) of 50° (Figure E}. This di- 
rection is approximately aligned to but slightly dif- 
ferent from the radio axis which has the PA of 75° 
(|Pentericci et al.ll2000D. Based on th e integral field spec- 
troscopy. iVillar-Martm et all (|2007l) found that the ex- 
tent of the Lya emission is ^9.0" x 7.5", and is misaligned 
by ^30±5° relatively to the radio structures. Therefore, 
the direction of the stretch of Ha emission is more or 
less similar to that of Lya emission. However, the size of 
Ha emission is ~ 4.5" x 3.0", and is much smaller than 
the Lya emission. It is possible that this smaller size of 
Ha emission is due to the fact that surface brightness 
of the outer region is less than the detection limit. To 
check the possibility, we compare the distributi on of Lya 
flux density given in IVillar-Martm et al.l (|2007l ) with the 
NB2315-if s image and find that the Ha emission can be 
detected at more than 2a if Lya/Ha ratio is constant 
over the whole region where Lya emission is detected. 
This thus means the ratio of Lya/Ha is not constant 
but larger in the outer region of the radio galaxy. Such 
smaller size of Ha emission compared to Lya emission is 
also seen in the r a dio galaxy MRC 1138-262 at z=2.16 
(jKurk et al.ll2002ft . iKurk et al.l (|2002f ) conclude that such 
fairly extended Lya emission is due to scattered radia- 
tion. Since we see a similar nature in the USS 1558-003 
radio galaxy, the existence of abundant scattering mate- 
rial around the radio galaxy trapped in deep potential 
wells and large Lya emission probably due to resonant 
scattering seem to be a common characteristic of HzRGs. 

Furthermore, we find that the radio galaxy also has ex- 
tended emission in B and r'-band, where there seems to 
be the contribution of Lya and He n emissions to some 
extend in the B and r'-band, respectively. However, we 
do not see any extended emission in the if s -band. Such 
extended emission in the rest-frame UV wavelengths, and 
the compactness of old stellar component may suggest 
that star forming activity is also currently on-going in 
the extended regions. In that case, such extended star 
formation is likely triggered by some physical processes 
associated to the radio activity. The jet-in duced star for- 
mation is one of the plau sible options fe.g. JBicknell et al] 
\2lM IHatchli^l2(nl . 

Many studies have been conducted to investigate the 
physical origin of nebular emission around the radio 
galaxy. Possible mechanisms include the photoioniza- 
tion by AGN, young stars and X-ray emission from 
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shock heated gas, and collisional ionizatio n by jet- 
induced shock (e.g- JMilev fc De Breuck 2008, and refer- 
ences therein). Although the ionization mechanism has 
not yet been completely understood, recent studies sug- 
gest that AGN could be the do minant source to ionize 
the surrounding nebular g as (e.g.. lHumphrev et al.ll2008l : 
iMilev fc De Breuckll2008l) . However, for radio galaxies at 
z ~ l. lBest et all (2000) found that the cause of nebular 
emission is dependent on the size of radio morphology 
based on the diagnostic with line ratios of nebular emis- 
sion. Small sources with a radio structure smaller than 
< 150 kpc tend to have emission originated from shock 
ionization, while larger sources prefer the mechanism of 
photoionization by AGNs. Since the radio observation 
of USS 1558-003 shows that the size of the radio galaxy 
is 9.2" (74kpc) (|Pentericci et al.ll2000fl . the shock ion- 
ization would be preferred. It is therefore interesting to 
investigate the origin of Ha emission of this radio galaxy 
more in detail. 

For this purpose, a diagnostic with ne bular emission in 
rest-fr ame optical wavelength is useful. iHumphrev et all 
( 2008) conducted optical and NIR spectroscopies where a 
slit was placed on the radio galaxy along the radio axis, 
and investigated the ratios of detected emission lines. 
They detected a lot of emission lines in the rest-frame 
wavelength range of A =1216-6585 A, and found that the 
radio galaxy has a broad Ha emission line with a ve- 
locity width of 12000 km s _1 , and that the amount of 
dust is negligibly small (A v ~ 0). Consequently, they 
concluded that the line ratios were best explained by 
AGN photoionization. However, they also find that the 
sources with high A v do not show evidence for jet-gas 
interactions, while the sources undergoing strong jet-gas 
interactions have low A v . The small amount of dust in 
this radio galaxy implies that the surrounding gas is un- 
dergoing interaction with the jet-induced shock. If this 
is indeed the case, the result is consistent with the pro- 
posed relation between the ionization mechanism and the 
radio size. It is suggestive therefore that the shock in- 
duced star formation is occurring to some extent as well 
as AGN photoionization, and thus both mechanisms are 
probably contributing to the large extent of the Ha emis- 
sion. It is difficult, however, to quantify the relative con- 
tribution of the two mechanisms, and it is beyond the 
scope of this paper. Integral field spectroscopy at NIR is 
essential to resolve the two-dimensional structure of the 
emission line regions and to fully understand the physical 
origins of nebular emission. 

6. CONCLUSIONS 

We have conducted a panoramic narrow-band imag- 
ing of Ha emitters (HAEs) in the proto-cluster candi- 
date around the radio galaxy USS 1558-003 at z = 2.53 
using NB2315 filter (A c = 2.313/im, AA = 0.027^m) in- 
stalled in MOIRCS on Subaru Telescope. This target is 
known as an over-dense region where distant red galax- 
ies (DRGs) are clustered (jKodama et al.l l2007). We have 
confirmed that this is indeed a rich proto-cluster in mak- 
ing with lots of star forming galaxies (HAEs) associated 
to the radio galaxy. We have mapped out the 2-D struc- 
ture of the proto-cluster, and investigated the star form- 
ing activities and the stellar mass content of this forming 
cluster. The main results we have found are summarized 
below. 
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Fig. 11. — Similar to Figure [5] but the size of the circles is now 
scaled with SFR, in the sense that larger symbols show higher 
SFRs of the Ha emitters. The red Ha emitters with J — Ks > 1.38 
(DRGs criterion) are also specified with red circles. The red filled 
dots show DRGs. 

(i) The proto-cluster is mainly composed of three con- 
spicuous groups of galaxies. One of them is sur- 
rounding the radio galaxy, and another is about 1.5 
Mpc (physical scale) away from the radio galaxy to 
the south-west, and the other is in between the two 
clumps. These groups show significant excess in the 
number densities of both HAEs and DRGs. Their 
close separations suggest that they would merge to- 
gether in the near future and grow to a single, more 
massive galaxy cluster at later times. 

(ii) A large fraction of the Ha emitters in this proto- 
cluster have SFRs higher than 100 M Q yr -1 , in- 
dicating that at z ~ 2.5, the progenitors of cluster 
early-type galaxies are vigorously forming in the bi- 
ased high density regions. Star formation activity is 
high everywhere irrespective of environment within 
the proto-cluster region, and the properties of indi- 
vidual HAEs show little environmental dependence, 
except that the HAEs in the densest clump may 
have slightly higher star formation rates compared 
to those in other regions. 

(iii) Most of the Ha emitters have blue colors, but some 
emitters have very red colors comparable to DRG 
(i.e., J — Ks > 1.38). Those red emitters are lo- 
cated on the fainter side of the red sequence on 
the color-magnitude diagram, except for the radio 
galaxy itself. Moreover, the red Ha emitters tend to 
be clustered in the three highest density clumps in 
contrast to lower- z clusters where similar red emit- 
ters are avoiding the cluster cores and preferentially 
located in the medium density regions or the out- 
skirts of the clusters. Since the red emitters are 
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likely to be dusty starburst galaxies in the tran- 
sitional phase, this result may indicate that some 
environmental effects, such as galaxy-galaxy inter- 
action, are at work on galaxies in the dense proto- 
cluster core at z — 2.53 and they are just changing 
their properties rapidly. 

(iv) The radio galaxy shows a large Ha halo extended 
over ^4.5" (i.e., 36 kpc), indicating that some ion- 
ization mechanisms (AGN, young stars, and shocks) 
are at work in the surrounding material around the 
radio galaxy. Such spatial extent of Ha emission is 
still much smaller than the Lya halo. This means 
that the Lya emission is severely extended by res- 
onant scattering. 

These results we present in this paper are all intriguing 
and there is no doubt that proto-clusters at z > 2 have 
high star forming activity and are in vigorously evolving 
phase. However, we need to investigate more samples 
of proto-clusters at z > 2 in detail in order to answer 
the question whether our results represent the universal 
properties of proto-clusters at z > 2 or showing the spe- 



cific characteristics of the USS1558 proto-cluster. Our 
on-going MAHALO-Subaru project has been surveying 
other proto-clusters as well as un-biasedly selected fields 
at similar redshifts, and enabling us to reveal the univer- 
sality or variation of properties among proto-clusters at 
z > 2. Such discussion will be presented in our future 
papers (e.g., Koyama et al. in prep.). 
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